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Abstract. We discuss a technique for mapping the syn- 
chrotron turnover frequency distribution using nearly si- 
multaneous, multi-frequency VLBI observations. The lim- 
itations of the technique arising from limited spatial 
sampling and frequency coverage are investigated. The er- 
rors caused by uneven spatial sampling of typical multi- 
frequency VLBA datasets are estimated through numeri- 
cal simulations, and are shown to be of the order of 10%, 
for pixels with the deconvolution SNR ~ 7. The fitted 
spectral parameters are corrected for the errors due to lim- 
ited frequency coverage of VLBI data. First results from 
mapping the turnover frequency distribution in 3C 345 are 
presented. 
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1. Introduction 



Information obtained with Very Long Baseline Interfer- 
ometry (VLBI) about radio spectra of parsec-scale jets 
and their evolution can be crucial for distinguishing be- 
tween various jet models. However, there are several as- 
pects of VLBI which impede spectral studies of parsec- 
scale regions. The reliability of spectral information ex- 
tracted from VLBI data depends on many factors includ- 
ing sampling functions at different frequencies, alignment 
of the images, calibration and self-calibration errors, a 
narrow range of observing frequencies, and source variabil- 
ity. The influences of all these factors must be understood 
and, if possible, corrected for, in order to reconstruct the 
spectral properties of parsec-scale jets consistently. 

Radio emission from the parsec-scale jets is commonly 
described by the synchrotron radiation from a relativistic 
plasma (e.g. Pacholczyk 1970). The corresponding spec- 
tral shape, S(v) oc v a , is characterized by the location 
of spectral maximum (S m , i^ m ) also called the turnover 
point, and by the two spectral indices, athick (for frequen- 
cies v -C I'm) and a t hin (for v > f m ). 



In many kiloparsec-scale objects, spectral index distri- 
butions have been mapped, using observations made with 
scaled arrays. In such observations, the antenna configu- 
rations are selected at each frequency in a specific way 
such that the spatial samplings of the resulting interfero- 
metric measurements are identical at all frequencies used 
for the observations. It is virtually impossible to use the 
scaled array technique for VLBI observations of parsec- 
scale jets made at different frequencies. The uneven spa- 
tial samplings of VLBI data at different frequencies result 
in differences of the corresponding synthesized beams, and 
can ultimately lead to confusion and spurious features ap- 
pearing in spectral index maps. 

In spectral index maps, the only available kind of in- 
formation is the spectral slope between the two frequen- 
cies. While sufficient for many purposes, this information 
can be misleading in the situation when the frequency of 
thespectral maximum lies between the frequencies used 
for spectral index mapping. In the ranges of frequencies 
between 1.4 and 43 GHz, frequently used for VLBI ob- 
servations, such a situation can be quite common. Using 
observations at three or more frequencies, it is possible 
to estimate the shape of the synchrotron spectrum, and 
derive the turnover frequency (frequency of spectral maxi- 
mum). Information about the turnover frequency can help 
to avoid the confusion which is likely to occur in spectral 
index maps. The turnover frequency is sensitive to changes 
of physical conditions in the jet such as velocity, particle 
density, and magnetic field strength. This makes it an ex- 
cellent tool for probing the physics of the jet in more detail 
than is allowed by analysis of the flux and spectral index 
properties of the jet. 

In this paper, we present a technique suitable for de- 
termining the turnover frequency distribution from multi- 
frequency VLBI data, and investigate its limitations and 
ranges of applicability. We discuss the advantages of us- 
ing the Very Long Baseline ArrayQ (VLBA) for spectral 
imaging. A general approach to imaging of VLBA data 
from nearly simultaneous, snapshot -type observations at 
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different frequencies is outlined in section 0. The effects of 
limited sampling and uneven itv-coverages are discussed 
in section ||. We provide analytical estimates of the sensi- 
tivity decrease, and use numerical simulations to evaluate 
the effect the uneven spatial samplings have on the out- 
come of a comparison of VLBI images at different frequen- 
cies. Alignment of VLBI images is reviewed in section |]. 
A method used for spectral fitting and determining the 
turnover frequency is described in section ||. Spectral fit- 
ting in th e ca se of limited frequency coverage is discussed 
in section 5J). The first results from the turnover frequency 
mapping are presented in section ^. 



2. Spectral imaging using VLBA data 

Many of the usual technical problems in making spectral 
index maps from VLBI data can be avoided by using the 
VLBA (see Zensus, Diamond, & Napier 1995). With the 
VLBA, it is possible to achieve array homogeneity, have an 
improved flux calibration, and make the time separation 
between observations at different frequencies negligible. 
The major remaining problems are image alignment and 
uneven spatial samplings of VLBI data taken at differ- 
ent frequencies. At an observing frequency v bs, the spa- 
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tial sampling of a baseline B — \/B\ 
by two antennas whose positions differ by Bx, By, Bz 
is characterised by the spatial frequency (e.g. Thompson, 
Moran & Swenson 1986) 



C = (f bs-B/c)sin6» . 



(1) 



where 9 is the angle between the baseline vector and the 
direction to the observed object. The quantity £ is often 
represented by its coordinate components u and v mea- 
sured in the spatial frequency plane (wi>-plane). The spa- 
tial sampling of a VLB array is described by the distri- 
bution of spatial frequencies accumulated during the ob- 
servation at all available baselines (itt>-coverage). In VLBI 
data at different frequencies, these distributions can differ 
significantly. 

To overcome, or at least reduce, the negative effect of 
uneven u ^coverages, the following scheme of observation 
and data reduction can be used: 

1) Quasi-simultaneous multi-frequency observations. 
The reliability of interleaved frequency observations can 
be demonstrated by Figure [j] which compares a VLBI im- 
age of 3C345 obtained from a full-track observation and 
an image made from simulated data that were sampled 
so that they represent the ut>-coverage achieved in a 5/15 
duty cycle observation (corresponding to an observation 
at three frequencies, with 5 minute -long scans at each 
frequency). Despite some loss of dynamic range and sen- 
sitivity, the extended structure is still well detected in the 
simulated image. The jet appearance in the simulated im- 
age is consistent with the jet structures seen in the actual 
map. 
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Fig. 1. Maps of 3C345 at 5 GHz (Zensus et al., 
in prep.) obtained from a real full-track global 
VLBI observation (top) and from a simulated 
interleaved-frequency VLBA+VLA observation (bot- 
tom). The contour levels are the same in both images: 
(-1,1,2,4,8,16,32,64,128,256,512) x 3.2 mJy. 



2) Careful choice of wavelengths. The choice of wave- 
lengths must be a compromise between the possibility 
of detecting the source structure and the VLBA sensi- 
tivity. Typically, at frequencies higher than 22 GHz, the 
requirements on brightness temperature limit the sensi- 
tivity. Also, there is a stronger dependence of the high 
frequency data on atmospheric instabilities. At frequen- 
cies lower than 2.3 GHz, many sources will become too 
complicated to warrant a successful structure detection 
without a full-track observation (see Table |l|) . 

3) Applying the phase-cal information for aligning the 
relative phases in separate frequency bands (Cotton 1995). 
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4) Improved amplitude calibration, due to frequent 
system temperature measurements and a relatively weak 
elevation dependence of the power gains of the VLBA an- 
tennas (Moran & Dhawan 1995). 

5) Applying appropriate wi>-tapering, in order to pro- 
vide matching wt>-ranges for data taken at different fre- 
quencies. 

6) Convolving data at different frequencies with the 
same, artificially circular beam. 

7) Using SNR and flux threshold cutoffs, in order to 
leave out the remaining sidelobe-induced artifacts. 

8) If strong sidelobe effects remain present, match- 
ing the wi>-coverages within certain wt>-ranges or over the 
whole wt>-plane. 

3. Spatial sampling 

Two effects are specific to the observing scheme described 
above: the decrease in sensitivity due to reduced uv- 
sampling at each observing frequency, and uneven spatial 
samplings at different observing frequencies. In this sec- 
tion, we investigate the effect of these factors on VLBI 
data at different frequencies. 

With a larger number of antennas, the overlapping 
parts of the ww-coverages at different observing frequen- 
cies constitute an increasingly larger fraction of the joint 
wi>-population. This results in decreasing the flux density 
level at which the confusion effects dominate the results 
of spectral index calculations. The confusion level can be 
lowered further by applying specific weighting schemes to 
the data (wy-weighting), and by matching the ranges of 
the ui>-coverages at both frequencies. In the spectral in- 
dex maps produced using the above procedures, confusion 
occurs at the flux level of about 0.3-0.5% of the mean 
peak flux density in the corresponding total intensity maps 
(Lobanov 1996). 

3.1. Longest sampling intervals and largest detectable 
structures 

For a given time sampling interval At, an estimate of the 
largest angular size, f2 max , of structures that can be de- 
tected on a baseline B can be calculated from the time- 
average smearing (Bridle & Schwab 1989). The sensitiv- 
ity reduction is greatest when the apparent motion of the 
source is perpendicular to the fringes associated with the 
selected baseline, and so we can assume a polar source 
and an East-West oriented baseline, in order to provide 
the most conservative estimates of f2 m ax- This gives 

max " w B cAt(B% + B*)V* ' { ) 

with u) e denoting the Earth angular rotation speed. For a 
feature located at the sky coordinates I, m with respect to 
the phase-tracking center (Thompson, Moran & Swenson 



1986), the corresponding average reduction in amplitude 
over a 12-hour period is 

< R At >=< I/I >« 1 - ^ {I 2 + m 2 sin 2 S) , (3) 

where S is the source declination. 

For every combination of wavelength and structure 
size, the left panel of Table gives the corresponding 
maximum allowed sampling interval [in minutes] between 
individual scans. Dots indicate that a structure remains 
unresolved. Italics highlight the 50% decrease of sensitiv- 
ity. The right panel of Table |l] gives the expected maxi- 
mum size of detectable structure [in mas], for all combi- 
nations of wavelengths and sampling intervals. The cal- 
culations have been done for the longest available VLBA 
baseline (B^ = 8600 km). Here we postulated Bz = and 
B 2 ^ = B 2 X + By, to provide more restrictive estimates. 

It follows from Table [l] that multi-frequency obser- 
vations can provide satisfactory structure and flux sen- 
sitivities for bright sources with intermediate (w 10 mas) 
extension. For such sources, full-scale spectral index map- 
ping with 5 minute-long scans can be done at frequencies 
lower than 43 GHz (0.7 cm). 

3.2. Simulations of multi-freguency VLBA data 

To study the effect of uneven vy-coverages on spectral 
imaging, we simulate visibility data at all frequencies 
available at the VLBA, using the routine "FAKE" from 
CIT VLBI package (Pearson 1991). At all frequencies, the 
simulated data are produced from the same "CLEAN" 
(Cornwell & Braun 1989) ^-component model of the struc- 
ture in the top panel of Figure [l[ To improve short-spacing 
coverage, we include one VLAnantenna in the simulations. 
The simulated bandwidth, W = 64 MHz, corresponds to 
one of the standard VLBA observing modes (128 Mb s -1 
data rate with 1 bit sampling; Romney 1992). We model 
the antenna efficiencies, rji (Crane & Napier 1989), us- 
ing the antenna sensitivities, Ki (Crane & Napier 1989), 
and the mean VLBA values, t/vlba and ifvLBA, given in 
Napier (1995). Then, for each VLBA antenna, the result- 
ing efficiency is: 

Vi = Wlba(^/^vlba) • (4) 

The system equivalent flux densities, SEFD, are cal- 
culated from the antenna sensitivities and system temper- 
atures, T sys : SEFDi = T aya<i /Ki (Walker 1995; Crane & 
Napier 1989). 

To make the simulated data as realistic as possible, we 
introduce four types of errors: the Gaussian additive noise, 
Ctherm, Gaussian multiplicative noise <r m , gain scaling er- 
rors again, and random station-dependent phase errors. 
Both <7 m and <7 ga i n are chosen to be at a 2% level, which 

2 The Very Large Array is operated by the National Radio 
Astronomy Observatory. 
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Table 1. Maximum sampling intervals [min] and maximum detectable structure sizes [mas] for a 8600 km-long baseline. 
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Table 2. Parameters of a typical VLB A and VLA antennas 
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is a good approximation of typical VLBA gain calibration 
errors. 

For a bandwidth of W [MHz] and integration time of 
T"int[s], the additive Gaussian noise can be calculated for 
each antenna, using the antenna zenith system tempera- 
tures, T sys , and antenna efficiencies r\. We have 

cthcrm = 5T sys / (rje p tDl nt ^Ti n tW) , (5) 

where D ant — 25 m is the antenna diameter, and e p t is 
the pointing efficiency. The pointing efficiency can be es- 
timated from the ratio of pointing errors, <7 p t, to the half- 
power beamwidth of an antenna at a given frequency. The 
typical non-systematic pointing errors of VLBA anten- 
nas are within 8-14" (Romney 1992), which results in 
€ pt ~ 85-99% at most of the VLBA observing frequen- 
cies. The overall parameters used in the data simulations 
are summarized in Table |^ for a typical VLBA (Wrobel 
1997) and VLA antennas. 

Using the parameters from Table ^, and adding the 
errors described above, we simulate VLBA visibility 
datasets that would be obtained in an observation with 
a 5/15 duty cycle corresponding to observing at 3 fre- 
quencies, with 5 minute long scans at each frequency. We 



f i t, 



Fig. 2. Baseline visibility amplitudes on the baseline be- 
tween VLBA antennas at Los Alamos and North Liberty. 
Top panel shows data from a real observation; simulated 
data are shown in the bottom panel. 



simulate a 17 hour-long observation of 3C345, with 30 
second averaging time for individual data points — similar 
to the typical duration and averaging time of real VLBI 
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Fig. 3. Fractional errors due to different Mi^coverages. 
The errors are plotted against the pixel flux scaled to the 
average self-calibration noise, a = 5 mJy. 

observations. The simulated data at 5 GHz are compared 
in Figure || with the data from a real VLBI observation, 
for a VLBA baseline Los Alamos - North Liberty. One 
can see that the noise levels are comparable in the real 
and simulated data. 

3. 3. Spatial sampling at different frequencies 

To study the effects of uneven wt>-coverages in VLBI data, 
we image the simulated datasets, following the procedure 
described in section p[ An image obtained from the simu- 
lated data at 5 GHz is shown in the lower panel of Figure |l|. 
Flux density and spectral index errors due to differences 
in spatial samplings can be estimated from comparison 
of the images made from the simulated data at different 
frequencies. In the ideal case, the flux ratio measured be- 
tween any two images should remain unity in every pixel, 
and the corresponding spectral index should be zero across 
the entire image. Since all images are produced from the 
same source model, we ascribe all deviations from zero 
spectral index to the errors due to different spatial sampl- 
ings and random errors, and choose to present these errors 
as a function of pixel SNR measured with respect to the 
self-calibration noise. The latter is taken to be equal to 
the largest negative pixel in the map, and is approximately 
5-10 times bigger than the formal RMS noise of the image. 
In the simulated data, the average self-calibration noise is 
« 5 m Jy. 

Figure || shows the fractional errors as a function of 
pixel SNR. We plot the results from all image pairs to- 
gether (5, 8, 15, 22, and 43 GHz data are used). The curved 
lines represent power law fits to the individual image pairs; 
the frequency ratio of each pair is given in the legend. The 
errors increase significantly at SNR< 7. We find that pixel 
SNR is the main factor determining the derived errors, al- 
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Fig. 4. Errors in spectral index due to different uv- 
coverages. The errors are plotted against the pixel flux 
scaled to the average self-calibration noise, a = 5 mJy. 

though the errors in pixels with comparable SNR tend to 
increase slightly at larger distances from the phase center. 
This increase however is considerably weaker compared to 
the increase of errors due to lower pixel SNR. 

The spectral index errors are shown in Figure ^ for the 
same image pairs. Similarly to the fractional errors, the 
magnitude of spectral index errors increases rapidly at low 
SNR. The main difference is that the errors become pro- 
gressively smaller at larger frequency separations, which 
follows obviously from the definition of the spectral index. 

From the error distributions shown in Figures we 
conclude that multi-frequency VLBA observations with 
the time sampling interval of 10 minutes at each frequency 
can be compared with each other, for pixels which are 
located at moderate (~ 10-15 mas) distances from the 
phase-tracking center, and have a sufficiently high SNR 
(> 5). Within these limits (and with the applied observing 
and data reduction strategy), the fractional errors should 
not exceed ^10%, a precision that can be sufficient for 
several purposes including spectral index and turnover 
frequency mapping in the nuclear regions of parsec-scale 
jets. Most of the large amplitude errors occur at jet edges 
where the effects of uneven spatial samplings are most 
pronounced. This effect is readily confirmed by Figure || 
in which the distribution of fractional errors is plotted for 
the 5-15 GHz image pair. The contours outline areas in 
which the errors are larger than 5%. These areas are con- 
centrated at the jet edges, and cover a fairly small fraction 
of the entire source structure. 

4. Image alignment 

Position measurements with the precision required for the 
alignment of VLBI images cannot be made without exten- 
sive absolute or relative astrometry observations. In quasi- 
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Fig. 5. Distribution of fractional errors in a 5-15 GHz im- 
age pair. Contours outline the regions with errors larger 
than 5%. 

simultaneous multi-frequency VLBI observations, the 



phase-referencing technique (Beasley fc Conway 1995) 
can be sufficient for the purpose of image alignment. If nei- 
ther of the aforementioned techniques is available, VLBI 
images are usually aligned by means of the position of 
compact core of the source. The core is likely to be lo- 
cated in an optically thick environment, and its position 
must depend on the observing frequency, v. According 
to Konigl (1981), the core is observed at the separation 
r CO re ^ v~ x l kr from the true jet origin. The term k r is close 
to unity (Marcaide et al. 1985, Lobanov 1997), so that 
r corc depends almost inversely on the frequency. There- 
fore, if aligned by the position of the core, VLBI images 
may contain systematic position offsets undermining spec- 
tral imaging. 

The frequency dependent shift of the core position can 
be deduced from comparison of observations made at close 
epochs, assuming that the superluminal features observed 
in the jet are optically thin and therefore should have their 
positions unchanged. In this case, the offsets between the 
component locations measured at different frequencies will 
reflect the respective shift of the observed position of the 
optically thick core. This approach has been successfully 
used in several studies (e.g. Biretta et al. 1986; Zensus 
et al. 1995; Lobanov 1997), and we deem it to be suffi- 
cient for the purpose of aligning VLBI images at different 
frequencies. 

5. Spectral fitting 

The frequency range covered by VLBI observations is very 
narrow (roughly, from 1 to 100 GHz, with most of the ob- 
servations done between 5 and 22 GHz). The turnover fre- 
quency, ^ m , of the synchrotron spectrum often lies outside 
the range of observing frequencies (see the sketch in Figure 
||). Because of the limited frequency coverage, a straight- 
forward application of the synchrotron spectral form to 
fitting may be ill-constrained (such as in the case B in 
Figure ^). To provide an estimate of v m in such cases, 
we first attempt to achieve the best fit of spectral data 
by polynomial functions, and then analyse local curvature 



true v nl 

Case B / 



fitted v„, 



■r. 
U 

Q 

x 



VLBI data 

X 



' fitted V™ 

/ 

/ 

/' Case A 



1GHz 



43GHz 



v \ 



-1-lOMHz 



Frequency 



-1000GHz 



Fig. 6. Sketch illustrating the problems existing in deter- 
mining the turnover frequency from VLBI data with in- 
sufficient frequency coverage. A homogeneous synchrotron 
source with isotropic pitch angle distribution is assumed. 
When the true turnover frequency lies outside of the range 
of observing frequencies, ad hoc information about low 
and high-frequency spectral shape is required, as well as 
the second order corrections based on the measured cur- 
vature of the observed part of the spectrum (see text for 
details). 

of the obtained fit. From the measured curvature, an es- 
timate (often only an upper limit) of can be found. 



Details of spectral fitting are given in section 5.4; appli 



cation of the local curvature for determining the turnover 



frequency is discussed in sections 5.5-5.6 



5.1. Basic synchrotron spectrum 

For an extensive coverage of synchrotron emission and 
its properties, we refer to the works by Ginzburg and 
Syrovatskii (1969), Pacholczyk (1970), and Ternov and 
Mikhailin (1986). In our calculations, we consider syn- 
chrotron emission from a homogeneous plasma with 
isotropic pitch angle distribution and power law energy 
distribution 71(7)^7 = n~ f0 ^/~~ s d'y, for electron Lorentz fac- 
tors 7l < 7 < 7h- In this case, it would suffice to describe 
the emission within the range of frequencies <C v <C ; 
where ^l,h are the low-frequency and high-frequency cut- 
offs given by: 



7l.h 



(6) 



where Q e is the electron gyro-frequency. Then, for a 
plasma with electron self-absorption, the spectral distri- 
bution of emission is (Pacholczyk 1970): 



1 — exp 



(7) 



where v\ is the frequency at which the optical depth, r s = 
1, and a t , a Q are the spectral indices of the optically thick 
and optically thin parts of the spectrum (with spectral 
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index denned by S oc v a ) . It is clear from equation [?] that 
at frequencies v <C v\. l v oc (v/vi) at ; and at frequencies 
v 3> V\\ J„ oc 1 — exp[— (vi/v) at ~ a °]. For a plasma with 
a homogeneous synchrotron spectrum at — 2.5. We will 
use the above at and the description given by (0) in our 
calculations. 



5.2. Fitting algorithm 

The main steps of spectral fitting can be summarized as 
follows: 

1) Make an approximate estimate of cutoffs in the spec- 
trum, on the basis of all available spectral information. A 
fairly good guess can be attained by taking the average 
measured turnover frequency and optically thin spectral 
index in the compact source. Based on these values, cal- 
culate the high and low frequencies at which the corre- 
sponding flux density is at an arbitrarily low level (we 
use >S cu tofT = 0.1 mJy). Add these spectral points to the 
measured data in each pixel, in order to ensure a negative 
curvature of the fitted curves, as required by the theoret- 
ical spectral form (R). 

2) Fit the combined pixel spectra by polynomial func- 
tions, allowing for limited variations of the cutoff frequen- 
cies, and aiming at achieving the best fit to the measured 
data points. From the fits, determine the basic spectral 
parameters: the turnover frequency, turnover flux density, 
and integrated flux (with integration limited to the range 
of observing frequencies). Estimate the errors from Monte 
Carlo simulations, using the distribution of the \ 2 param- 
eters of the fits to the simulated datasets. 

3. Calculate the local curvature of the fitted spectra 
within the range of observing frequencies. Compare the 
derived curvature with the values obtained from analyti- 
cal or numerical calculations of the synchrotron spectrum. 
Derive the corrected value of the turnover frequency, by 
equating the fitted and the theoretical spectral curvatures. 

4. Fit the data with the synchrotron spectral form 
described by (]?]), and using the corrected value of the 
turnover frequency. 

The above procedure has been applied to spectral data 
obtained from modelling VLBI images of 3C 345 by ellip- 
tical Gaussian components, and combining the models at 
different frequencies (Lobanov & Zensus 1998). 

5.3. Spectral cutoffs 

Because the transition between high-frequency (y > vf) 
and low-frequency (y < vf) spectral regimes determined 
by equation ^ is very sharp, the spectrum is determined 
by the synchrotron self-absorption at frequencies v <C v\ , 
and by the electron energy distribution at v ^> v\. Esti- 
mates of the typical turnover frequency and energy spec- 
tral index in the jet of 3C345 based on our own calcula- 
tions and on the results from Rabaca and Zensus (1994) 
give vi 10 GHz and s « 2.6. Using these values, we esti- 



mate the low-frequency, ~ 1 MHz, and high-frequency, 
~ 1000 GHz, cutoffs in the spectrum, and use these val- 
ues for the spectral fitting. However, the cutoffs cannot be 
well defined, and may change during the evolution of the 
jet emission. To account for this effect, we allow 15% vari- 
ations of the cutoff frequencies so as to achieve the best 
fit to the data. 

5.4. Spectral fits 

In order to calculate the spectral parameters of jet emis- 
sion, we combine the component fluxes measured at fre- 
quencies y±, ...,z/Ni and add the cutoff information. The 
resulting spectral dataset S s , is represented by the flux 
densities, 

{S s } = {S(^l),<SVi),---,SVn),SVh)}, 
and their respective variances, 

Ws} = {°sK), ctI(^i), •-, cts(^n), ctsOh)} , 

where S'(^'l) = S(vn) = 0.1 mJy is the cutoff flux density 
level. 

By varying the spectral dataset, we then produce sim- 
ulated spectral datasets 



Ssym — V(S s )\{al(u)} 



(8) 



assuming the Gaussian distribution of the errors in flux 
density measurements 



<f>(S) 



1 



27r<7„ 



(9) 



where and cr„ are the measured flux density and its 
variance. 

For each simulated spectral dataset, we apply the 
M th -order (M < N — 1) polynomial fitting by the ba- 
sis functions X k {v) = v k , and minimize the \ 2 parameter 
of the fit 



N 1 
E ^2 



M 



(S, 



sym J i 



k=0 



(10) 



to obtain the least squares fit to the data, and derive the 
polynomial coefficients 



M 
fe=0 

Here aj.k and f3k are given by 

JV 

CLj,k 



(11) 



i=l 
N 



Xj(vi)X k (vi) 
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Synchrotron Spectrum Calculation 
Fitted/True Turnover Frequency Rati 





Fig. 7. Theoretical curvature, k, of the homogeneous syn- 
chrotron spectrum with spectral index a. The £ axis de- 
notes the ratio of the frequency at which the curvature is 
calculated to the turnover frequency. 

From the fit to the spectral dataset, we derive the basic 
parameters of the synchrotron spectrum: the integrated 
flux, 



S\nt — 



M 

Aj dv . 



(12) 



the turnover frequency, v m , and the turnover flux density, 
S m : 



dv 



=> S m ,v„ 



(13) 



We then analyse the distribution of the % 2 parameters 
from the fits to all simulated datasets. From this analysis, 
standard deviations at the 3er confidence level are calcu- 
lated for S-mt, S m , and u m . 

5. 5. Curvature of the fits 

The local curvature of spectral fits is: 



d 2 S 
dv 2 



1 



AS 

d7 



-1/3 



(14) 



The mathematical details of calculations are summarized 
in Appendix. If the spectral index, a , and the local cur- 
vature, k , of a polynomial fit are determined at a fre- 
quency v Ql then the turnover frequency, v m , can be es- 
timated from the adopted theoretical synchrotron spec- 
trum. Using the derived a and k , we determine the ra- 
tio £(a ,K ) = v m /v Q , from the adopted spectral form 
described by (Q). The corresponding turnover frequency is 
then 



Levs= 1.00DOE+00'( 0.500, 0.700, 0.850, 
1 .000, 1 .200, 1 .500, 2.000, 3.000, 5.000, 
10.00, 20.00, 30.00, 50.00, 75.00) 



Fig. 8. Curvature correction coefficients derived from sim- 
ulating synchrotron spectra and subsequent fitting them 
by the polynomial functions described in section 5.4. 

Figure Q relates the curvature k of the homogeneous syn- 
chrotron spectrum described by (Q) to the ratio £,(a , n ). 
For frequencies increasingly deviating from the turnover 
frequency (for which £ = 1), the curvature, n, becomes 
progressively smaller, thereby limiting the ranges of ap- 
plicability of the corrections described by (|l5|). For data 
covering the frequencies from v\ to z/n, we expect the cor- 
rections to give reliable estimates for the turnover frequen- 
cies lying within the 0.05^! < v m < range. 

5. 6. Numerical estimates of the curvature corrections 

An alternative method to estimate £(a ,K ) is to deter- 
mine it numerically, by simulating the synchrotron spec- 
trum with given turnover frequency and spectral index, 
and fitting it by the polynomial functions used in section 
Km We have performed such calculations, using the spec- 



tral form described in section 5.1, and covering a range of 



turnover frequencies and spectral indices. The results are 
shown in Figure || in which the ratio of fitted values to 
the theoretical values of the turnover frequency is plotted 
against the synchrotron spectral index. The contours show 
different values of £. One can see that the determined ra- 
tios do not depend strongly on the spectral index. Figure 
H can be used for the same purpose of correcting the fit- 
ted turnover freque ncies obtained through the procedure 
described in section 5.4. 



v m = v ^(a , k ) . 



(15) 



6. Turnover frequency mapping 

On the basis of the method described above, we have de- 
veloped a code for mapping the turnover frequency distri- 
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bution from multi-frequency VLBA data. Fitting is per- 
formed in every valid pixel of the image; the validation 
is based on clipping the pixels with low flux density or 
low SNR. Flux density errors are estimated from the noise 
level and flux density gradients in the total intensity maps. 
For each pixel, we average the values of pixels within se- 
lected bin, and add, in quadratures, the averaging stan- 
dard deviations to the estimated noise level. This results in 
slightly increased errors for pixels in the areas with steep 
flux density gradients, providing more conservative error 
estimates. The use of the gradients for error estimation is 
optional, and can be turned off by setting the bin size to 
1. 



The output of the mapping procedure can be the 
turnover frequency distribution, turnover flux density dis- 
tribution, integrated flux distribution, or total intensity 
map at a given frequency. The last option allows us to pre- 
dict, from the fitted spectral shape, the expected source 
structure at any frequency within the range of observing 
frequencies of the maps used for the spectral fitting. This 
can also be used for testing the quality of the spectral fit, 
by comparing the predicted and observed images at the 
same frequency (given that the observed image was not 
used for producing the above spectral fit). 



6.1. Mapping the turnover frequency distribution in 
3C345 

The blazar 3C345 (z = 0.594, Hewitt & Burbidge 1993) 
is a strongly variable core-jet type source with a compact 
core responsible for most of the source radio emission, and 
a curved, parsec-scale jet (Zensus et al. 1995) containing 
enhanced emission regions (bright components) travelling 
along curved trajectories, with speeds of up to 20 c (Zensus 
et al. 1995). Synchrotron spectra of the core and the near- 
est bright components are often peaked around 10 GHz 
(Lobanov & Zensus 1998), and show a remarkable evo- 
lution. The emission from the core and the components 
is believed to be produced by condensations of highly- 
relativistic electron-positron plasma injected in the jet, 
and losing their energy first through the inverse-Compton 
mechanism (Kellermann & Paulini-Toth 1969; Unwin et 
al. 1997), and later on due to the synchrotron emission 
from adiabatically expanding relativistic shocks (Wardle 
et al. 1994; Zensus et al. 1995). 

The turnover frequency procedure was applied to the 
multi-frequency VLBA observation of 3C 345 made on 
June 24, 1995. The source was observed at 5, 8.4, 15.4, 
and 22.2 GHz. At each frequency, there was roughly one 
5 minute scan made every 20 minutes. After the corre- 
lation, the data were fringe-fitted and mapped in AIPS^ 

3 Astronomical Image Processing Software developed and 
maintained by the National Radio Astronomy Observatory 
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Table 3. Parameters of the VLB A maps 



1 


2 


3 


4 


5 


6 


7 


8 


^obs 


Stot 


•Speak 




ocst 
^noisc 


Beam 


u«-range 


w-taper 


[GHz] 


[Jy] 


[Jy/bm] 


[Jy/bm] 


[mjy/bm] 




[MA] 


[MA] 


22.2 


7.210 


4.138 


-0.015 


4.4±0.6 


0.75x0.63, 9. 3 


0-574 


150 


15.4 


7.321 


4.439 


-0.015 


3.9±0.4 


0.84x0.68, 5.°5 


0-440 


150 


8.4 


7.607 


4.173 


-0.007 


2.4±0.2 


1.16x0.87, -20.°3 


0-240 


150 


5.0 


7.103 


3.803 


-0.018 


4.1±0.9 


1.57x1.20, -4.°1 


0-150 


150 



Notes: 1 - observing frequency; 2 - total CLEAN flux; 3 - peak flux density; 4 - minimum flux density; 5 - measured noise; 6 
- major axis, minor axis, and position angle of the tapered beam; 7 - uu-range of the data; 8 - half-power taper size. 



4x10 3 2x10 3 -2x10 3 -4x10' 




Relative R.A. [arcsec] 

Fig. 10. Turnover frequency distribution in the extended jet of 3C345. The central region is saturated, for better 
representation of the turnover frequency variations in the jet. The contours are drawn at 0.1, 0.15, 0.2, 0.3, 0.7, 1, 2, 
5, 8, 10, 12, and 15 GHz. All values below 5 GHz should be regarded as upper limits. 



and DIFMAP (Shepherd 1993). The data were tapered 
at 150 MA, and the maps were produced with a circular 
restoring beam of 1.2 mas in diameter. The core shift with 
respect to the reference frequency (22.2 GHz) was applied 
to the data at 5, 8.4, and 15.4 GHz. The magnitude of 
the core shift for the data at 15.4 GHz was determined 
from the fit 7- CO rc <x l/ -i-04±o.i6 j w ] lereas f or the data at 
5 and 8.4 GHz the measured values were used (Lobanov 
1998). The resulting maps are shown in Figure B; the main 
characteristics of the maps are given in Table~p|. Marked 
in the maps are the source core "D" and jet component 



"C7" which dominated the source emission at the epoch 
of observation. 



The turnover frequency map produced from the VLBA 
maps shown in Figure |^ is presented in Figure [l^. Figure 
pM| shows a map at 11 GHz obtained from the spectral fit 
to the combined VLBA data at 4 frequencies. One can see 
that the main features in the predicted 11 GHz image are 
consistent with the structures seen in the original VLBA 
maps. 
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2 - 3C345 11 GHz VLBA 24/06/95 




Relative R A. (milliaresec) 



Fig. 11. Map of 3C 345 at 11 GHz obtained from the spec- 
tral fit. The restoring beam and contour levels are the 
same as in Figure 10. Spectral profiles in Figure 13 are 
taken along the horizontal line crossing the nuclear region 
of the source. 

6.2. Nuclear region 

In Figure [w], there are two regions of higher turnover 
frequency in the nucleus of 3C 345 oriented nearly trans- 
versely to the direction of the jet. These regions match 
the locations of the core and C7 fairly well. The increased 
turnover frequency may indicate that the emission is com- 
ing from a shocked plasma. The transverse extension is 
then consistent with strong shocks that are likely to be 
oriented almost perpendicularly to the jet direction. Fig- 
ure ^| shows spectral profiles made along the horizontal 
line crossing the center of the core (horizontal line in Fig- 
ure 0. The core and C7 are both visible in the turnover 
frequency profile. The turnover flux distribution is very 
smooth and peaks almost precisely at the center of the 
core. From the turnover frequency and turnover flux dis- 
tributions, we can derive the profile of magnetic field in 
the central region using the relation (Cawthorne, 1991) 

B(r) = C ^r 4 5 m 2 , (16) 

where Co is the proportionality coefficient. Co can be de- 
termined empirically from the estimates of the absolute 
position (r co 



5 pc) and magnetic field (B c 



0.3 G) 



of the core at 22.2 GHz (Lobanov 1998): 

Co = ^corc'S'm.corc^'m^corc ~ 1-2 ' 10 ° , (17) 

for the measured 

^m^orc — 5.6 Jy and k'nijCorc — 

15.1GHz. 

Equation [l(] expresses the magnetic field strength due to 
the compression that the plasma has undergone during 
shock formation. Therefore, the magnetic field also de- 
pends on the strength of the underlying magnetic field in 
the location of the jet where the shock is formed. We pos- 
tulate that the underlying magnetic field i? a mb oc r~ m , 
and consider the cases, with m = 1 and m = 2. The jet 
is assumed to have a constant opening angle (f> = 2.4 
(Lobanov 1998). For the magnetic field in an arbitrary 
pixel p, formula [lq yields 



20 



15 



10 







turnover frequency [GHz] 






turnover flux density [Jy] 






magnetic field [G] 






\ m=1 






^ 






m=2 ( \ 




/ \ / 


\. \ ' -A 
ri Vi\ 









Br, 



2 p( r p/ r core) 4 m [G] 



(18) 



-5 -4 -3 -2 -1 
Relative R.A. [mas] 

Fig. 12. Profiles of turnover frequency, u m , turnover flux, 
S m , and magnetic field, B{r), along the line AS = — 3mas 
crossing the center of the core (horizontal line in Figure 
12). The underlying magnetic field decreases along the jet 
as r~ m . 



In this formula, is measured in GHz, S m is in Jy, and 
r is in parsecs. The resulting magnetic field profiles are 
plotted in Figure [l^. The magnetic field rises sharply, close 
to the outer edge of C7. This can signify the amount of 
plasma compression in the shock. The increased magnetic 
field on the opposite side (particularly visible in the B oc 
r~ 2 profile) may reflect a larger electron plasma density 
near the jet origin. In the relativistic jets, the case B oc r _1 
is expected to be more likely. A somewhat high value of 
the magnetic field in C7 (-B™7 X ~ 8.5 G) in this case may 
also be caused by possible errors in the estimates of the 
core magnetic field. However, the derived shape of the 
magnetic field profile is consistent with C7 being a strong 
shock embedded in the jet of 3C345. 

6.3. Extended jet 

Almost everywhere in the extended jet shown in Fig- 
ure [l0[ the turnover frequency is lower than 5 GHz, pos- 
ing a problem for both the spectral fitting and assessing 
the results from the fits — we therefore resort to regard- 
ing all values of v m < 5 GHz as upper limits. Appar- 
ently, there are no strong shocks dominating the extended 
jet of 3C345, or their turnover points may have evolved 
rapidly due to strong adiabatic cooling. Because the de- 
rived turnover frequencies are too low, we cannot make 
quantitative statements about the physical conditions in 
the extended jet. Observations at lower frequencies (1.6, 

1.4, 0.6, 0.3 GHz) are required for a better understanding 
of the turnover frequency changes in these regions. With 
the available data, we can only make general comments 
about the gradients observed in the turnover frequency 
map. The bright patterns elongated along the jet ridge line 



12 



A. P. Lobanov: Spectral imaging with VLBI data 



may indicate the presence of an ultra-relativistic channel 
inside the jet (e.g. Sol et al 1989). The extended patterns 
seen in the jet at oblique angles to the ridge line resemble 
the patterns of Kelvin- Hclmholtz instabilities (see Hardee 
et al. 1995, for the results from 3D simulations of the 
KH-instability driven jets). As has been noted above, the 
turnover frequency is exceptionally sensitive to the varia- 
tions of plasma speed and density. Therefore, the observed 
patterns may reflect the velocity gradients and/or den- 
sity gradients existing in the jet perturbed by the Kelvin- 
Hclmholtz instability. However, the low frequency data 
are needed for making a better substantiated conclusion 
about the observed gradients. 

7. Summary 

In this paper, we have covered several methodological and 
scientific aspects of studying synchrotron spectrum of the 
parsec-scale regions in AGN. The main conclusions can 
be stated as follows: 

1) We have discussed a technique that can be used 
for mapping the turnover frequency distribution and ob- 
taining spectral information from multi-frequency VLBA 
data. A feasibility study shows that multi-frequency 
VLBA observations can be used for spectral imaging and 
continuous spectral fitting. 

2) Multi-frequency VLBA observations made with up 
to 10 minute separations between the scans at each fre- 
quency can provide a satisfactory spatial sampling and 
image sensitivity for sufficiently bright sources with inter- 
mediate (~ 10-15 mas) structures. The fractional errors 
from comparing the data at different frequencies should 
not exceed 10% for emission with SNR> 7, in this case. 

4) A procedure for broadband synchrotron spectrum 
fitting has been introduced for mapping the distribution 
of spectral parameters of radio emission from parsec-scale 
jets. Corrections based on the local curvature of the fit- 
ted spectra are introduced, in order to compensate for 
the incomplete frequency coverage in cases where the true 
turnover frequency is outside of the range of observing 
frequencies. 

5) From a 4-frequency VLBA observation of 3C345, 
the first map of the turnover frequency distribution are 
produced. The maps indicate possible locations of the rela- 
tivistic channel and strong shock fronts inside the jet. The 
magnetic field distribution derived from the turnover fre- 
quency and flux distributions is consistent with the plane 
shocks existing in the immediate vicinity of the source 
core. The extended emission appears to have a very low 
turnover frequency for which the existing data do not war- 
rant a good estimate, limiting the conclusions to deduc- 
ing certain information from the gradients of the turnover 
frequency which are visible in the extended jet. The ob- 
served gradients are consistent with the patterns of veloc- 
ity distribution and density gradients typical for Kelvin- 
Hclmholtz instabilities propagating in a relativistic jet. 



A more detailed study, with observations made at lower 
frequencies, is required for making conclusive statements 
about the nature of the observed gradients of the turnover 
frequency. 
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Appendix: Calculation of the local curvature of the 
fitted and theoretical spectral forms 

Linear spectral forms 

With the fitted polynomial coefficients ao,...a3, we can 
write a linear form of the fit as: 

S(v) =C exp(r), (Al) 
with Co = exp(ao) and 

3 

t = ^(^(lnz/)* . 

And the derivatives used in ( |l4] ) are given by the following 
formulae: 



dS dr 
— = C exp(r) — 
dv dv 



d 2 S 
dv 2 



C 



, ,d 2 r dr 

ex P r X~2 + ex P 2r X" 
dv £ dv 



dr 1 

— = — \a\ + 2a 2 In j/ + 3a 3 (\nv) 2 ] 
dv v 



(A2) 
(A3) 
(A4) 



d 2 T i 

-prr = ^[3a 3 (ln^) 2 + (2a2-6a 3 )ln^ + (ai-2a2)](A5) 

d^ z v^ 



The power-law fit is given by: 



S{v) = Ci 



v\ 



1 — exp 



with the corresponding derivatives: 



dS_ 
dv 



= C l 



a t ^/" i (a t -l) t 

v *— V a t v y ' f v 

dv 



d^ " ClW d^ 



dv 

+ (a t - l)a t f v ] 



(A6) 



(A7) 



(A8) 
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/„ = 1 - exp 



A = a Q — a t 



exp 



d^ = ^^P-l)(vj exp 
Logarithmic spectral forms 

Following the same considerations, the logarithmic form 
for the polynomial fit and its derivatives are: 



2A-2 



(A9) 
(A10) 
(All) 



dS'iog 



i=0 



d(lni/) 



= a\ + 2a 2 v + Sa^v 



d 2 S 



log 



= 2d2 + 603^ 



d(ln^) 2 
And for the power-law fit: 
Slog(v) = InCi + a t (lnv - lni/i) + In/,, 



T7i \ = A — exp 

d(lni/) 



d(ln^) 2 \v! 

dg'iog 1 dfv 
d(lni/) ~ /„d(lni/; 



1- - 











exp 









+ «t 



d 2 S ; 



log 



d(lni/)2 / 2 



d 2 /„ / d/„ 



d(ln^) 2 Vd(ln^) 
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